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ABSTRACT 

We present time-resolved spectrophotometry of the pulsating DA white dwarf G 29-38. 
As in previous broad-band photometry, the light curve shows the presence of a large 
number of periodicities. Many of these are combination frequencies, i.e., periodicities 
occurring at frequencies that are sums or differences of frequencies of stronger, real 
modes. We identify at least six real modes, and at least five combination frequencies. 
We measure line-of-sight velocities for our spectra and detect periodic variations at the 
frequencies of five of the six real modes, with amplitudes of up to 5kms _1 . We argue 
that these variations reflect the horizontal surface motion associated with the g-mode 
pulsations. No velocity signals are detected at any of the combination frequencies, 
confirming that the flux variations at these frequencies do not reflect physical pul- 
sation, but rather mixing of frequencies due to a non-linear transformation in the 
outer layers of the star. We discuss the amplitude ratios and phase differences found 
for the velocity and light variations, as well as those found for the real modes and 
their combination frequencies, both in a model-independent way and in the context of 
models based on the convective-driving mechanism. In a companion paper, we use the 
wavelength dependence of the amplitudes of the modes to infer their spherical degree. 

Key words: stars: individual (G 29-38) - stars: oscillations - white dwarfs 



1 INTRODUCTION 

Pulsating white dwarfs show brightness variations at a range 
of frequencies, with amplitudes of up to a few per cent. They 
appear to be normal white dwarfs, except for the fact that 
they lie in a special temperature range, or instability strip, 
in which they are unstable to pulsations. There are instabil- 
ity strips for all three of the main white dwarf types, DA, 
DB and DO (atmosphere dominated by H I, He I and He II, 
respectively), around 12 000, 23 000 and 100 000 K, respec- 
tively. 

Pulsating white dwarfs offer the prospect of physical 
insight via a number of different methods. One is asteroseis- 
mology, the determination of a white dwarf's internal struc- 
ture through the comparison of observed periodicities with 
model predictions of stellar eigen-modes. Another method 
analyses which modes are present, what their amplitudes 
are and how the amplitudes vary. The aim is to learn about 
the physics of the upper layers of the star, by determining 
the way in which modes interact with different parts of the 
star and with each other. 



Robinson, Kepler & Nather (1982) showed that white 



dwarfs pulsate in non-radial gravity modes, and that the 
observed brightness variations result from perturbations of 
their temperature, not of their geometry or gravity. They 
also estimated the expected variations in the line-of-sight 
velocity and the line profile, and concluded that these would 
be too small to measure with then available instrumentation. 
As we will show below, this statement no longer holds. 

The pulsations were originally thought to be excited 
in an ionisation zone, as in classical pulsators such as the 
Cepheids. The dynamics in DA and DB pulsating white 
dwarfs are rather unlike those in classical pulsators, however, 
in that for the white dwarfs the surface convection zone can 
respond almost instantaneously to gravity-mode pulsations, 
because the turnover times are much s horter than the mo de 
periods. In a series of papers, Brickhill fll983| , |l990| , |l991a[ |b|) 
has described the response in detail, and concludes it leads 
to significant "convective driving", which will be sufficient 
to overcome the damping in the radiative interior for modes 
with periods of order of or shorter than the thermal response 
time of the convection zone. Brickhill's partly numeri cal re- 
suits have recent ly bee n confirmed analytically (Wu 
Goldreich & Wu |l999a| ). 
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The DA and DB pulsating white dwarfs often show 
power at frequencies that are the sums or differences of fre- 
quencies with stronger power. Based on the density of the 
theoretical gravity-mode spectra in the relevant frequency 
range, it has been argued that these "combination frequen- 
cies" do not represent real eigen-modes, but are instead pro- 
duced by nonlinear mixing of pulsation signals of real eigen- 
modes ("parent modes") in the outer layer of the star. Some 
support for this idea comes from observed correlations in 
power. 



Brickhill ( 1992a , b) has suggested that the mixing is a 
by-product of the interaction between the pulsations and the 
convection zone. Because of the relatively slow thermal re- 
sponse of the convection zone, flux variations take some time 
to appear at the surface, and they are reduced in amplitude. 
But the thickness of the convection zone is modulated by the 
pulsation, so that the signal is not only diminished and de- 
layed, but also distorted. It is the latter distortion which 
translates into the power seen at combination frequencies. 

Both to be able to make an asteroseismological analysis 
and to verify the pulsation models, it is essential to identify 
periodicities with eigen-modes, i.e., obtain their radial or- 
der n, spherical degree I and azimuthal order m. For some 
DB and DO variables, it has been possible to obtain rea- 
sonably sound identifications from a direct comparison of 
observed frequencies with model predictions. However, for 
many other stars, including all DA variables, there are ei- 
ther too few modes, or the mode structure is too variable. 



Robinson et al. (1995) implemented a method to iden- 
tify spherical degree i, which combined broad-band measure- 
ments of mode amplitudes at optical and UV wavelengths. 
It makes use of the wavelength dependence of limb dark- 
ening, which causes the pulsation amplitudes to vary with 
wavelength in a manner that depends on I, Based on its 
success, we were motivated to implement a variation of this 
method which uses high signal-to-noise time-resolved spec- 
trophotometry at optical wavelengths only. As will become 
clear, our initial data on the star G 29-38 showed variations 
not only in flux but also in line-of-sight velocity. The broad- 
ened scope, as well as the increased complexity in modeling 
the flux variation as a function of wavelength, made us de- 
cide to present the amplitudes and phases of the modes as 
a f uncti on of wavelength in a separate paperj^j (Clemens et 



1998 , hereafter Paper II). 



al. t 

In this paper, we focus on the implications of the line- 
of-sight velocity variations. We first briefly review the basic 
properties of G 29-38 (Section ^), and then describe the ob- 
servation and data reduction in detail (Section ^). In Sec- 
tion ^, we analyse the periodicities in the light curve, and in 
Section ^ we search for corresponding variations in the line- 
of-sight velocity curve. We discuss the implications of our 
results in Section |(| both in a model-independent way and 
in the context of the convective driving models described 
above. 



We are far from having exhausted the information in our data; 
we will gladly make it available in digital form to anyone inter- 
ested in analyzing it further. 



2 G 29-38 

G 29-38 (ZZ Psc; WD 2326+02 9) wa s discovered as a vari- 
able by Shulov & Kopatskaya (11974) and studied in some 
detail by McGraw and Robinson (1975). It is bright (V ~ 

11800 K, 

Bergeron et al. 



13.0), at the c ool en d of the instability strip (T ■■ 

cf. Paper II) and shows large-amplitude 



I!)!).") 



modulations (up to a few per cent). The sour ce rem ained rel- 
atively obscure until Zuckerman & Becklin (1987) reported 



an infrared excess, which is still unex plained (for recent de- 
velopmen ts, see Kleinman et al. 



1994; Koester, Provencal 



Shipman 1997 ). In 1987, it was first subjected to sc rutin y 
with the Whole Earth Telescope (WET; Winget et al. h990h. 



Unlike for the variable DBs and DOs, however, the observa- 
tions - while of excellent quality - yielded more questions 
than answers. Partly, this is because of the complex and 
variable frequency spectrum of G 29-38: typically, dozens of 
periodicities are seen, with many occurring at frequencies 
that are differences or sums of the frequencies of other peri- 
odicities. 



Kleinman et al. (1995) discussed an extensive collec- 
tion of data sets, obtained over many seasons. They find 
that many periodicities are variable and unstable from year 
to year, but some do recur. By considering the whole en- 
semble, they could identify with some certainty which pe- 
riodicities were associated with real modes and which with 
combination frequencies. They found nineteen real modes, 
the majority of which very likely have 1 = 1. However, due 
to the vagaries of the mode characteristics - especially the 
enigmatic, confusing behavior of multiplets, showing differ- 
ent and changing splitting - their results are suggestive, but 
not definitive. 



3 OBSERVATIONS 

Spectra of G 29-38 were taken on 1996 November 19 usin 
the Low Resolution Imaging Spectrometer (Oke et al 
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on the 10-m Keck II telescope. From 4:42 till 9:25 UT, a 
continuous series of 700 12 s exposures was made. Only 100 
pixels in the spatial direction (corresponding to 21.5 arcsec), 
binned by a factor 2, were read out through two amplifiers, 
leading to a read-out and re-set time in between exposures 
of about 12 s. The 600 line mm -1 grating was used, covering 
the range 3450-5950 A at 1.2 A pix . The weather started 
photometric, and we used a 8.7 arcsec slit to make use of 
that. Hence, the seeing of 1.2 arcsec determined the wave- 
length resolution of ~7A. In principle, it would have been 
best to keep the slit at the parallactic angle in order for re- 
fraction not to influence the wavelength scale, but this would 
have required non-standard guiding and we decided to keep 
the slit at position angle 0° instead. 

At the end of the series, a bank of cloudlets came in, 
causing the last 55 spectra, as well as one spectrum some 
ten minutes earlier, to have a reduced source and increased 
sky rate (the moon was up, so these are easy to identify). 
We have not used any of these spectra in our analysis. 

The series of exposures was followed by fifteen 4 s in- 
tegrations on the spectrophotometric flux standard G191- 
B2B, using the same setup. Each of the series was fol- 
lowed by Hg/Kr wavelength calibration and halogen flat- 
field frames. 
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The reduction of all spectra was done as follows. The 
frames were bias-corrected separately for the two amplifiers 
using the overscan regions, and then corrected for the gain 
difference between the two amplifiers as determined from 
halogen spectra. The sky background was estimated by fit- 
ting parabolae for each position in the dispersion direction, 
excluding 21 binned pixels (~ 9") centred on the stellar spec- 
trum. The latter pixels were used to estimate the stellar flux 
using an optimal weighting scheme similar to that of Home 
( 1986 ). Note that we did not flat-field the data. The reason 
is that we had been unable, due to an instrument failure 
shortly after taking the exposures on G191-B2B, to obtain 
a large enough number of halogen spectra to construct a 
flat-field that matched the signal-to-noise ratio of our ob- 
servations, especially in the blue. Since the flat fields that 
we did have appeared very smooth, we decided to forego 
flat-fielding altogether. The resulting spectra look smooth 
(Fig. Q), except for a few pixels that are easy to identify. 
Note the presence of metals - Ca II A3933 a nd M g II A4481 
- uncovered only recently by Koester et al. (1997). 

Accurate flux calibration was not possible, since the 
G 191-B2B spectra were taken while the bank of cloudlets 
was still passing. We selected the three most exposed spec- 
tra of the series of fifteen acquired, as these have compara- 
ble count rates and hence may not be influenced too much, 
making the relative flux calibration more reliable^. From 
these three, we derived the instrumental response usi ng th e 
model fluxes provided by Bohlin, Colina & Finley ( [l995[ ), 
and the at mosph eric extinction curve of Beland, Boulade 
& Davidge (1988). Typical examples of relatively high and 
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low-flux spectra, as well as the average spectrum, are shown 
in Fig. §. 

Special care was taken to correct the wavelength scale 
for instrumental effects due to flexure and atmospheric re- 
fraction. Both effects can be seen in the pixel positions on 
the detector as a function of time of the O I A5577 sky line 
and the Balmer lines (Fig. ^a). Here, the positions were de- 
rived by cross-correlation of the spectra with the average of 
those spectra taken with the slit position angle less than 20° 
away from the parallactic angle. The break at t — i ~ 2000 s 
in the Balmer-line positions is due to guiding being stopped 
and re-started (to take a guider image). 

The flexure is described well by a cubic fit to the sky-line 
positions. After correction for the flexure, as well as for the 
jump due to switching the guider off and on, the positions 
of the Balmer lines vary linearly with what is expected from 
refraction (sin a tan £, where a is the difference between par- 
allactic and position angle, and £ the zenith distance). The 
different slopes indicate that our guiding wavelength was ap- 
proximately 7000 A, as expected for the CCD camera used. 
We used this value to calculate the wavelength corrections 
for the whole spectrum. 

Since our spectra were taken through a wide slit, small 
movements of the star will result in random apparent wave- 
length shifts for individual spectra. The effect can be esti- 
mated by measuring the shifts along the slit from the spatial 
profiles (averaged over some number of pixels along the dis- 



t Perhaps surprisingly, our estimated flux in the average spec- 
trum compares quite well with the obse rved V magnitude of 
13.05 ± 0.05 (22 ± 1 mjy; McCook & Sion |l987|). 



20 



10 



n — i — I — i — i — i — i — I — i — i — i — r 



l — I — i — i — i — r 




4000 4500 5000 

x(A) 



5500 



Figure 1. Average of all spectra and examples of spectra taken 
at times of high and low flux (times —7444 and —7080 s in 
Fig. ^b). Apart from the Balmer lines (marked), weak lines due 
to Ca II A3933 and Mg II A4481 are also present (see insets; size 
40 A by 2mjy). The other irregularities in the average spectrum 
are due to pixel-to-pixel sensitivity variations (see text). 



persion direction). We find a root-mean-square variation in 
these of 0.084 binned pixels, corresponding to 0'.'036. Assum- 
ing the variations in the dispersion direction are similar, one 
thus expects variations of about 0.2 A, or 14kms _1 at H7. 

Apparent brightness variations can result from varying 
slit losses due to seeing-related changes in the width of the 
stellar images. From the spatial profiles, we find 10% root- 
mean-square variations around the average full width at half 
maximum of 3 binned pixels (1'.'2). Using images taken with 
LRIS for other projects, we estimate that the slit losses with 
an 8'.' 7 x 9'.'0 aperture are about 1% on average. Assuming 
that the slit losses scale approximately quadratically with 
the seeing, one thus expects apparent brightness variations 
of the order of 0.2%. This is smaller than the Poisson noise 
in individual wavelength elements, but dominates in the av- 
erage light curve discussed below. It is also much large r than 
variations due to scintillation (~0.012%; Young 1967). 



4 PERIODICITIES 

In Fig. ^b, we show the continuum light curve, constructed 
by averaging the 5200-5500 A region. The light curve shows 
the presence of a number of periodicities, and is similar to 
light curves observed before (see Section ^). The Fourier 
transform is shown in Fig. Most of the power is at fre- 
quencies below 6mHz. An idea of the background level can 
be obtained from the 10-20 mHz region (see inset in figure), 
where little evidence for any periodicity is seen. In this re- 
gion, we find a root-mean-square amplitude of about 0.025%. 
The corresponding root-mean-square variation is half that, 
and one infers root-mean-square variations in the fluxes of 
about i^/AdZT x 0.025% = 0.32% (where A" data = 644 is 
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Figure 2. (a) Variations as a function of time in the pixel po- 
sition on the detector of the O I A5577 sky line (bottom curve; 
offset by —5 pixels), H/3 (top curve at right) and He (middle curve 
at right). The smooth variations are due to flexure in the spectro- 
graph and differential refraction in the atmosphere, (b) Fractional 
brightness variation (top curve) in the line-free wavelength region 
5200— 5500 A, and residual variation (bottom curve, offset by 0.9) 
remaining after taking out the eleven periodic variations with 
amplitudes larger than 0.55% and a constant term (Fig. ^a; Ta- 
ble h). (c) Line-of-sight velocity variations determined from fits to 
the H/3, FTy and H<5 lines (top curve), and the residuals (bottom 
curve, offset by -100 km s -1 ) remaining after taking out varia- 
tion at the eleven periods present in the light curve, as well as a 
cubic term. 



the number of data points), more or less consistent with the 
estimate made above based on varying slit losses. We believe 
essentially all the signal seen in the main panel of Fig. |^a, 
which is well above 0.025%, reflects stellar variability. 

We identified the periodicities sequentially. We started 
by measuring the approximate frequency of the highest 
peak and fitting the data to a function of the form C + 
A cos(2tv ft' — 4>), where t' = t— 6:52:08 UT is the time relative 
to the mid-observing time, solving for pulsation frequency 




2000 4000 
F CuHz) 



6000 



Figure 3. Fourier transform of the light curve (a) and line-of- 
sight velocity curve (b) . The eleven strongest peaks in the Fourier 
transform of the light curve are indicated and labeled. The pluses 
indicate the amplitudes derived from fitting eleven cosine waves 
(plus a constant for the light curve and a cubic term for the radial- 
velocity curve) . The bottom curve in each panel shows the Fourier 
transform of the residuals, offset from zero for clarity. The insets 
show the Fourier transform at higher frequencies. 



/, phase <fi and amplitude A, plus a constant offset C. Next, 
we identify the highest remaining peak in the Fourier trans- 
form of the residuals and include another cosine wave in the 
fit, etc. For combination frequencies, we held the frequen- 
cies fixed at the corresponding combinations of parent-mode 
frequencies. The decision whether a peak corresponded to a 
combination frequency or to a real mod e, was generally made 
consulting the work of Kleinman (1995); in all cases, the am- 



plitude of the combination frequency was the smallest of the 
peaks involved. Continuing as long as the next largest peak 
had not been influenced greatly by removing the previous 
peaks, we found eleven periodicities with amplitudes larger 
than 0.55%. 

Of the eleven periodicities found, there appear to be 
six real and five combination frequencies. These are listed 
in Table |l| and the corresponding peaks are identified in 
Fig. C ompa ring our periodicities with those found by 
Kleinman ( |L995| ), we find that most periodicities have been 
seen before. The exceptions are: (i) F6: power at 430 s has 
only been seen as an apparent harmonic of a strong 860 s 
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Table 1. Strong Pulsation Frequencies in G 29-38 



Real mode 


P 


/ 


A L 




*L 




A v 


4> v 


Ry 


A$ v 




00 


OuHz) 


(96) 




(°) 




(kins" 1 ) 


(°) 




(°) 


Fl 


614 


1628.0 ±0.9 


3.17 ± 


08 


-168 ± 


2 


5.4 ± 0.8 


-116 ± 8 


17 ± 3 


52 ± 8 


F2 


818 


1223.0 ±0.9 


2.71 ± 


08 


132 ± 


2 


2.3 ±0.8 


176 ± 19 


11 ± 4 


44 ±19 


F3 


653 


1531.2 ± 1.8 


1.46 ± 


08 


36 ± 


3 


2.4 ±0.8 


59 ± 19 


17 ± 6 


23 ± 19 


F4 


776 


1288.7 ±3.9 


0.90 ± 


08 


40 ± 


6 


2.0 ±0.8 


50 ±22 


27 ±11 


10 ±23 


F5 


283 


3528.8 ±4.9 


0.59 ± 


08 


55 ± 


8 


1.6 ±0.8 


131 ± 27 


12 ± 6 


76 ±28 


F6 


430 


2327.6 ±4.9 


0.62 ± 


08 


-86 ± 


8 


0.4 ± 0.8 




4± 9 




Combination 


P 


/ 


A L 




*L 






4>v 


Rc 


A* c 


frequency 


00 


(/.Hz) 


(%) 




(°) 




(kms -1 ) 


(°) 




(°) 


F2 + F1 


351 


2851.1 


0.97 ± 


08 


-33 ± 


5 


0.6 ±0.8 




5.7 ±0.5 


3± 5 


F2 - Fl 


2469 


405.0 


0.67 ± 


08 


34 ± 


7 


1.1 ± 0.8 




3.9 ±0.5 


-26 ± 7 


2F2 


409 


2446.1 


0.95 ± 


08 


-131 ± 


5 


0.6 ±0.8 




12.9 ± 1.4 


-35 ± 7 


F3 + F1 


317 


3159.2 


0.58 ± 


08 


-129 ± 


8 


1.4 ±0.8 




6.2 ± 1.0 


3± 9 


F3 - Fl 


10323 


96.9 


0.66 ± 


09 


159 ± 


7 


1.6 ±1.1 




7.1 ± 1.0 


3± 8 



combination frequencies 



Notes. For the real modes, Ry = (Ay /27t/)/Al in units of 10kmrad _1 %~ 1 and A<l>v = *v — *L- F° r the 
/ = /* ± f 1 , with for difference frequencies i and j chosen such that / is positive, Rc = 
(with riij the number of possible permutations of i and j; unity for i = j, 2 otherwise) and A<E>c = 
± "l?^)- The uncertainties in the parameters for the light curve should be seen as indicating relative 
uncertainties only, as the light curve is not described well by the eleven modes fitted and the deviations from the fit 
are definitely not distributed normally (see text and Fig. Mb). The estimates given here were obtained by multiplying 



In 



the observational errors by such a factor that \ 



2 

rod 



1. For the velocity curve, the estimates of the uncertainties 



should be reliable. The unc ertainty distributions on amplitude and phase are not normal, but those on A cos <E> and 
Asin<E> are (see Section 5.1). 



mode; (ii) F3±F1: power at Fl (614 s) and F3 (653 s) has 
not been seen simultaneously before, and hence neither has 
this combination (although power at 317 s due to a differ- 
ent combination has been observed); and (iii) F2— Fl and 
F3— Fl: at 2469 and 10322 s, these periods are longer than 
any generally considered (the identification of F3— Fl should 
be considered tentative here too) . Absent is power at 400 s, 
which has been present in almost all observations so far. 

From the residuals, it is clear that the eleven modes 
do not account well for the observed light curve (with our 
estimate of 0.2% root-mean-square uncertainties, we find a 
formal Xicd — 50). We have tried continuing the identifi- 
cation process (see appendix), but fail to account for all 
the variability. Therefore, we decided to focus on the eleven 
modes for which we have full confidence. We note, however, 
that because of the additional power and the resulting non- 
normal deviations, the formal errors on the parameters are 
not very meaningful. To get an idea of relative uncertainties, 
we have scaled up the errors on the data points such that 
the reduced x 2 equals unity. 

A puzzling aspect of our results is that the three largest 
periodicities appear to have frequencies related to each other 
by simple integer ratios: fi : fi ~ 3 : 4, /2 : fa ~ 4 : 
5. Both are correct to within 0.2%. After identifying more 
modes and including these in the fit (see appendix), the 
ratios become more accurate (to within 0.1%). Furthermore, 
another mode, F8, is identified for which fs : fi — 2 : 3, 
again to within 0.1%. Similar close coinciden ces h ave not 
been found in the observations of Kleinman (1995) and it 



may be that they are in some way due to our relatively 
short time span, perhaps in combination with the badness 
of the fit, especially near the light maxima (see Fig. |b). 
On the other hand, these ratios may represent unexplored 
physics underlying the pulsation spectrum. For the present 
time, we will ignore possible complications arising from this. 



5 VELOCITY VARIATIONS 

The positions of the Balmer lines on the detector, as derived 
by the cross-correlations described in Section |§| appear to 
show modulation on time-scales similar to those shown by 
the light curve (see Fig. ^a). This modulation might reflect 
apparent variations in the line-of-sight velocity of G 29-38. 
In order to study this in more detail, we determined Doppler 
shifts for all the s pect ra by fitting the H/3, H7 and H<5 lines 
(details in Section 5.2 below). The resulting line-of-sight ve- 
locity curve is shown in Fig. ^3 (upper curve). The Fourier 
spectrum, shown in Fig. bp, shows one strong peak, at a 
frequency of 1634 ± 5/iHz, i.e., consistent with that of the 
strongest periodicity in the light curve. 

Inspired by the close correspondence between the fre- 
quencies of these signals, we looked for power at frequen- 
cies corresponding to the other ten periodicities in the light 
curve and found small peaks at some but not all of them 
(Fig. ^) . We determined their amplitudes and phases by fit- 
ting the velocity curve with a combination of eleven cosine 
waves with frequencies fixed at those determined from the 
light curve (plus a third-order polynomial to account for 
slow variations). The results are listed in Table 1 and indi- 
cated in Fig. ^b; the residuals are shown in Fig. 7c (bottom 
curve). Note that unlike for the light curve, the fit to the 
line-of-sight velocities is formally acceptable (x?cd — 1) an( A 
we expect our error estimates to be reliable. 

Could these line-of-sight velocity variations reflect pul- 
sational motion? The pulsations are g-modes and the associ- 
ated motion is largely horizontal (e.g., Dziembowski 1977). 
The observed line-of-sight velocity, i.e., the physical veloci- 
ties projected on the light of sight and integrated over the 
stellar disk, will thus be dominated by motion near the limb. 
The observed amplitude will be smaller than the physical ve- 
locity amplitude by an amount depending on the spherical 
degree £ and the direction of the pulsation axis relative to 
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the observer. Robinson et al. ( |1982| ) estimated the expected 
line-of-sight velocity amplitudes for typical pulsation ampli- 
tudes and found that they should be a few kms -1 . This is 
comparable to the amplitudes we find. 

If we suppose that the line-of-sight velocity variations 
reflect horizontal motion associated with the white dwarf 
pulsation, it is very interesting that all five apparently sig- 
nificant detections (> 2a) are for periodicities we identified 
as real modes above; no significant detections are found for 
the combination frequencies. Before we turn to the implica- 
tions, however, we need to verify that our error estimates are 
reasonable and the velocity detections significant; and if so, 
whether the inferred line-of-sight velocities reflect intrinsic 
velocity variations rather than some systematic effect in our 
measurement procedure. We will address these two points 
in turn. 

5.1 Significance of the detections 

Looking at Fig. ^b, all peaks but that associated with Fl 
look insignificant. Indeed, they would be formally insignifi- 
cant if one did not know the frequency at which to look in 
advance. However, the pulsation frequencies are determined 
accurately from the light curve. In the absen ce o f system- 
atic effects (considered separately in Section 5.2), a rough 
estimate for the significance of a peak with a certain am- 
plitude at a frequency specified in advance can be made by 
counting the number of other peaks in the velocity spec- 
trum at or above that amplitude. For instance, there are 13 
peaks with amplitudes above 2 kms -1 at frequencies below 
7mHz (Fig. ^b). With a frequency resolution of ~ 64/iHz, 
the total number of independent frequency elements in this 
frequency range is 109, and thus the probability of finding a 
> 2 kms -1 peak by chance at a given frequency is roughly 
1:8. One would estimate 1:14 if one counts peaks up to the 
Nyquist frequency (21 mHz; 23 peaks over 2 kms -1 in 322 
[iVdata/2] frequency elements). Thus, F2, F3 and F4, all of 
which have amplitudes in excess of 2 kms" 1 , are likely real. 
There is no other peak as large as Fl, hence the probability 
that it is a chance coincidence is < 1 : 322. 

In order to obtain more quantitative estimates, we ran 
Monte-Carlo simulations, in which we fitted scrambled ve- 
locity curves in the same way as we fitted the observations. 
The scrambled curves were created by first taking out the 
slow third-order variation in the measured Doppler shifts, 
and then re-distributing the resulting velocities randomly 
over the measurement times. By counting the number of 
artificial data sets for which the fit produces an amplitude 
larger than that observed, one can infer the likelihood of a 
false detection. The likelihood distributions for all periodici- 
ties but F3— Fl are consistent with each other (because of its 
long period, F3— Fl is covariant with the third-order pol 
normal, and its likelihood distribution is wider). In Fig. 
we show the cumulative distribution of the amplitudes for 
all modes but F3— Fl derived from one thousand simula- 
tions. We also indicate the observed amplitudes. We infer 
that Fl, F2, F3 and F4 are detected at the > 99.9, 98, 98 
and 95% confidence levels, respectively; F5 is marginally de- 
tected (~88% confidence); and the signals at all other peri- 
odicities are not significant. 

As mentioned above, the fit to the line-of-sight veloc- 
ities is formally acceptable (Xred — !)■ The errors, domi- 
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Figure 4. Cumulative distribution of line-of-sight velocity am- 
plitudes derived from simulations with scrambled data. The two 
dotted lines indicate the amplitudes at which a measurement at 
a frequency specified in advance has a 90 or 95% probability of 
being real (i.e., not due to chance). The crosses indicate the am- 
plitudes derived from the fit to the observed line-of-sight velocity 
curve (Fig. |^). One sees that the detections of real modes Fl, 
F2, F3 and F4 are significant (J; 95% confidence), that of F5 is 
marginal, and that of F6 is not significant. No significant signal is 
measured for any of the combination frequencies (F3— Fl is not 
shown, because its probability distribution is different; see text). 



nated by the star wandering around in the slit (Section y|), 
should be more or less normally distributed, and thus our er- 
ror estimates likely reasonable. Indeed, from the simulations 
one infers similar numbers: ga = [§(^4 2 )] 1/ ' 2 — 0.8 kms -1 
(where the factor | accounts for the random phase), and 
as = (180/7r)(0.8kms -1 /Av) deg. The exception is F3-F1, 
for which a a = 1.5 kms -1 , indicating that the uncertainty 
for this periodicity is underestimated in the least-squares fit 
(Table hi). This reflects the large covariance with the third- 
order polynomial. 

While the error estimates should be reliable, we caution 
that the uncertainty distributions are not expected to be 
normal. This is a consequence of our choice of parameters 
Ay and $v- For instance, from Table [l], the uncertainty on 
f&v is 28° for F5. The 2er uncertainty, however, is far larger 
than 56° : indeed, at 95% confidence $v can have any value, 
as the detection of velocity variations at F5 is significant 
at < 95% confidence. The uncertainties on Ay cos <E>v and 
Ay sin $v are normally distributed, and one can use these 
to determine reliable uncertainty estimates for Ay and $v 
at different confidence levels. Note that the same holds for 
Al and $l, but the more significant detections make this 
less of an issue. 



5.2 Possible systematic effects 

Above, we have shown that at frequencies at which the flux 
of the star varies strongly, significant variations are also 
present in the line-of-sight velocity curve. The question that 
remains is whether these variations reflect the intrinsic mo- 
tion in the white dwarf's outer layers expected from the es- 
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timates of Robinson et al. (1982), or whether they might be 



biased or perhaps even result from inadequacies in the way 
we determined the velocities. For instance, through changes 
in continuum slope, the measured Doppler shift might be 
influenced by the variations in temperature with pulsation 
phase. 

For any systematic effect related to the flux variations, 
the line-of-sight velocity variations have to be either in phase 
or in anti-phase with the light variations, i.e., the phase dif- 
ference between velocity and light A"l>v = $v — $l = 
or 180 deg. This is not what we observe: for Fl, A<I>v = 
52 ± 8 deg. Thus, no more than Ay cos A$v = 3.3 km s -1 
could be related to systematic effects of this kind. The com- 
ponent out of phase with the light variations has an ampli- 
tude of j4vsinA$v = 4.2 km s -1 , which we can be confi- 
dent is real given the results of the Monte-Carlo simulations 
(Fig. §). 

While the above shows that at least part of the line-of- 
sight velocity signals are intrinsic to the star, the component 
of the variation that is in phase or anti-phase with the light 
variations may still be influenced. In order to get a handle 
on the importance of possible systematics, we tried different 
methods to determine the line-of-sight velocities, looked at 
differences between the velocities inferred for different lines, 
and tried fitting artificial data generated with the models 
described in Paper II. 

We found that cross-correlation was very susceptible to 
systematics. For instance, the line-of-sight velocity ampli- 
tudes derived from the cross-correlation discussed in Sec- 
tion ^ (Fig. ^a) are systematically different for different 
Balmer lines, both in phase and in amplitude. Most likely, 
this reflects the influence of the varying temperature and 
hence continuum slope, which is not accounted for in the 
cross-correlation procedure. To get unbiased velocities, one 
needs to normalize all spectra, but it is not clear how to do 
this properly. 

We therefore focused on line-profile fitting, of H/3, H7 
and H<5 (in the wavelength intervals 4622-5102, 4220-4540 
and 4040-4200 A, respectively) . To each line, we fitted a sum 
of a Lorentzian and a Gaussian profile, with their central 
wavelengths forced to be the same, and a continuum repre- 
sented by either a linear or a parabolic function. We tried 
fitting to logF in addition to fitting to F directly, the ra- 
tionale being that since absorption is a multiplicative rather 
than an additive process, a sum of components can represent 
log F better than F. 

For all types of fits, the resulting line-of-sight veloc- 
ity curves were fitted with the combination of eleven cosine 
waves. The results were systematically different from those 
inferred from the cross-correlation velocities, with ampli- 
tudes generally being somewhat smaller, and phase differ- 
ences between light and velocity being further away from 
zero. The systematic differences between lines were smaller, 
but not absent. No apparent effects are present in the de- 
rived phases for the different periodicities, but the ampli- 
tudes increase systematically from H/3 to H7 to H<5. Also, 
the zero points (i.e., systemic velocities) derived for the dif- 
ferent lines are very different. This is probably unrelated, 
however, as we find similar differences in systemic velocity 
for our flux standard G 191B2B, which has much narrower 
lines. Instead, it likely reflects problems with the wavelength 
calibration, possibly related to the way the light of the arc 



lamp is projected on the slit (T. Bida, 1997, private com- 
munication). Such an effect should not affect the results for 
the periodic variations in velocity. 

This leaves the systematic differences in the amplitudes 
found for the individual Balmer lines. We can discard the 
possibility that these are intrinsic to the star based on a sim- 
ple analytic derivation: while the intrinsic velocity and light 
variations are not independent, their dependency is only a 
second-order effect and hence very small (Paper II). Most 
likely, therefore, the observed systematic effects in ampli- 
tude are still related to our way of fitting, perhaps simply 
because our choice of Lorentzian plus Gaussian is an inade- 
quate model for the observed line profiles. 

To get a better handle on the problem, we fitted a num- 
ber of artificial data sets, produced using the procedures 
described in Paper II, for different settings of temperature 
amplitude, horizontal surface velocity amplitude, phase dif- 
ference between light and velocity, and spherical degree^. 

For all artificial data sets and fitting methods, the re- 
sults for the different lines showed systematic differences of 
comparable magnitude to those seen in the real data. A set 
with zero input temperature variation gave the same veloc- 
ity amplitude and phase for all three lines, as should be the 
case, but different systemic velocity, with H7 being close to 
zero, and H/3 and H<5 being negative and positive, respec- 
tively, by an amount that depended on the fitting method. 
The amount was smaller for fitting log F than for fitting F 
directly, as expected (—13, —1, -l-llkms" 1 for H/3, H7, H<5, 
respectively). 

For data sets with non-zero temperature variations, H7 
again reproduced the input values best, in the sense that 
the systemic value was still closest to zero, and the velocity 
phase delay closest to the input value. H/3 and H<5 tended to 
lower and higher amplitude and to larger and smaller phase 
delay, respectively, with the values for both amplitude and 
phase being more or less symmetric around those for H7. 
From a set with no velocity variations put in, we find that 
for a temperature variation corresponding to 2.8% modula- 
tion at 5500 A, the apparent velocity shifts in H/3 and H<5 
are about 0.6 km s -1 in amplitude. The systematic differ- 
ences with H7 for sets including velocity variations are con- 
sistent with this. The amplitude should be compared with 
the observed radial velocity amplitude of ^Skms" 1 and the 
uncertainty of 0.8 kms -1 . 

Given the results on the artificial data and taking into 
account that the observed velocity signals have low signal- 
to-noise, we decided to use the average of the H/3, H7 and 
H<5 velocities, determined by fitting to log F (with a linear 
slope for the continuum). The resulting line-of-sight velocity 
curve is the one discussed above and shown in Fig. ^c. It has 
uncertainties dominated by the effects of the star wandering 
around in the slit. Given the simulations with the artificial 
data, we feel reasonably confident that systematic effects in 
the fitted amplitudes and phases will be smaller than the 
quoted uncertainties derived from the least-squares fit (as 
listed in Table [l]). For completeness, we note that the result 



• In principle, we could have tried fitting the models directly 
to the observations. We decided not to do this both because the 
models do not reproduce the observed variations all that well, and 
because it would make verifying our results much more difficult. 
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would be similar if one were to use the velocities derived 
from H7 only. 

5.3 Velocity summary 

We conclude that we have detected line-of-sight velocity 
variations which are due to intrinsic processes in the white 
dwarf, and which have amplitudes comparable with those 
expected from the integration over the visible disk of th e 
surface motion due to the pulsations (Robinson et al. 1982). 



We have shown that our estimates of the uncertainties are 
robust and that systematic effects, while problematic, most 
likely will be within these uncertainties. Out of the eleven 
strongest modes present in the light curve, five out of six of 
the real modes show significant modulation, while none of 
the combination frequencies does, even though F2+F1 and 
2F2 have light amplitudes similar to those of the weaker real 
modes F4 and F5. 



6 DISCUSSION 

In this section, we discuss the observational results for the 
real modes and the combination frequencies in some detail. 
For each group, we describe the general properties we re- 
gard as important, mention what can be understood from 
geometrical arguments alone, and discuss the results in the 
context of convective-driving models. 



6.1 Real modes 

We mentioned already that only the real modes show sig- 
nificant modulation in the line-of-sight velocities. As mea- 
sures of the relative amplitudes and phases of velocity and 
light, we use Ry = {Ay /2irf)/A L and A<E> V = "l>v — $l 
(see Table |l|). Here, Ay and Al are the observed line-of- 
sight velocity and fractional light amplitudes, <I>v and 
the corresponding phases, and / the frequency of the mode. 
The ratio Ry does not depen d on azimuthal order m and 



aspect (see Dziembowski 



1977 



Paper II). The rationale for 



including / is that with this definition Ry depends only 
on spherical degree I for adiabatic pulsations. For adiabatic 
pulsations, one also has A<E>v = 90°, i.e., velocity maximum 
arrives a quarter cycle after light maximum. 

From Table [j], we see that the values of both Ry and 
A&y show no significant variation among the different real 
modes. All values of A$v are between and 90°, i.e., ve- 
locity lags light, but by less than a quarter cycle. 

The largest value of Ry is found for mode F4. While 
the difference with the values for the other modes is not 
formally significant, we note that a larger value is consis- 
tent with this mode having 1 = 2 and the others 1=1, 
as inferred in Paper II from the wavelength dependences of 
the pulsation amplitudes. Compared to an £ = 1 mode, the 
light variations for an I = 2 mode are subject to stronger 
cancellation, while the surfac e velo cities are actually easier 
to observe (see Dziembowski 1977). 



In the interior of a white dwarf, pulsation associated 
with a real mode should be largely adiabatic; maximum 
flux coincides with maximum compression and is followed 
a quarter cycle later by matter moving away at maximum 
velocity. This will change as the mode propagates upward 



and enters the outer regions, where non-adiabatic effects be- 
come important. Using the observed relations between light 
and line-of-sight velocity variations, we can therefore hope 
to infer properties of the surface layers. 



In the c onvect ive-driving picture (Brickhill 1992a; Gol- 
dreich & Wu 1999a), the main source of change is the surface 



convection zone. As this region bottles up heat when it is 
heated, the flux perturbations at the photosphere of the star 
are diminished in magnitude and delayed in phase compared 
to those entering the bottom of the convection zone. In con- 
trast, the horizontal velocities associated with the pulsation 
are nearly independent of depth inside the convection zone , 
because of the strong turbulent viscosity (Brickhill 199C ; 
Goldreich & Wu 1999b). Thus, values of Ry should be larger 



than those expected for adiabatic pulsations, and values of 
A$ v smaller than 90° . 

The effects are expected to become more prominent 
with increasing values of lutc, where u) = 2-rvf is the an- 
gular frequency of the mode and tc is the so -called ther- 
mal adjust ment time (Goldreich & Wu 1999a; quantity D 
in Brickhill 1983), which is a few times the thermal time at 



the bottom of the convection zone. With increasing mode 
frequency, therefore, Ry should increase and A&y should 
deviate more and more from 90°, tending towards 0°. 

Indeed, the observed values of A<E>v are smaller than 
90° and those for Ry larger by a factor 2-3 than expected 
if one were to neglect non-adiabatic effects (e.g., i?v,ad — 6 
for mode Fl). Unfortunately, the observations do not allow 
one to verify the predicted trends, but it may be worth not- 
ing that, if anything, Af&v shows a trend opposite to that 
predicted, i.e., less of a deviation from 90° for shorter peri- 
ods. We remind the reader that the errors in A$v are not 



distributed normally (Section 5.1) 



It is difficult to make a more quantitative comparison, 
as this requires evaluating the thermal adjustment time tc, 
which depends on details of the convection. Ideally, one 
would like to use Ry and A«I>v to estimate tc and thus 
constrain the convection properties. Our data only allow a 
rough determination: a 800-s mode in a star wit h tc = 300 s 
has Ry ~ 16 and A$ v =d 50° (Wu & Goldreich |l998|) , sim- 
ilar to what is observed for mode F2. 

An independent constraint to tc can be inferred from 
the longest-period mode that is present. This is because to 
excite a gravity-mode, the convection zone needs to extend 
deep enough for its thermal response time to be of orde r of or 
longer than the mode p eriod, i.e., ujtc > 1 (Brickhill 



198J 



Goldreich & Wu 1999a). The longest-period mode is F2 at 



818 s (or F7 at 1106 s; see Appendix), and thus tc <; 130 s (<; 
180 s) , consistent with the above value. Consistent estimates 
of tc are also obtained using combination frequencies (see 
below). 

6.2 Combination frequencies 

The lack of line-of-sight velocity signals at the frequencies 
which we associate with combinations of real modes, con- 
firms the notion that combination frequencies do not reflect 
physical pulsation, but rather mixing of the parent-mode sig- 
nals by a non-linear transformation occurring in the outer 
layers of the white dwarf. 

We quantify the combination frequencies using an am- 
plitude ratio Rc = A\^ /mjA^A^ and a phase difference 
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Here, A^ 3 and <k^ 3 are the ob- 



served amplitude and phase of a combination at frequency 
f l ±/ J of real modes i and j (with frequencies / ! and f 3 , am- 
plitudes A]^ and A 3 ^, and phases $l an d THj is the num- 
ber of possible permutations: mj = 1 for i = j and riij = 2 
for i ^ j; i and j are chosen such that difference frequencies 
are positive. We include ny to ensure that the harmonic of 
a single mode and a combination of two modes with almost 
identic al properties would lead to the same Rc (cf. Brickhill 
1992b). Similar de finit ions are used for the 3-mode combi- 



nations (see Table Al); below, we focus on combinations of 
two modes. 

The observed values of Rc and A$c (Tables [j] and Al) 
vary significantly from mode to mode, although part of the 
variation is due to only a few outliers (especially 2F4, which 
we will discuss in more detail below). There appears to be 
some correlation between the value of Rc and the ampli- 
tudes of the parent modes, with Rc increasing for decreas- 
ing A^f . The phase differences A$c cluster around zero 
(within ~45°); this reflects the sharp maxima and shallow 
minima in the light curve. 

If the combination frequencies indeed result from non- 
linear mixing in the upper layers of the star, then inde- 
pendent of a specific theory one expects a combination fre- 
quency to have at any point on the surface a flux amplitude 
proportional to the product of its parent-mode amplitudes at 
that point. The surface distributions of the parent modes are 
described by spherical harmonics. Therefore, the surface dis- 
tribution of the combination frequency will be proportional 
to the product of spherical harmonics, which can in turn be 
described by a linear superposition of spherical h armo nics 
with different I (see, e.g., Abramowitz & Stegun 
adopt their notation). 

In the following discussion, we will assume all parent 
modes have only m = components. We consider the more 
general case later. First, we study the first harmonic of a 
single parent mode. The first harmonic of a parent mode 
with (£,m) = (1,0) has a surface distribution described 
by YiYi, which can be decomposed into a sum of Y ° and 
y 2 ° distributions. The first harmonic of a parent mode with 
(£, m) = (2, 0) has a surface distribution that is the sum of 
three components, Y Q ° , Y§ and Y® . The numerical factors 
in front of the Y Q ° terms are the same for both decompo- 
sitions. When integrated over the visible hemisphere, the 
(2,0) parent mode will suffer stronger cancellation than the 
(1,0) parent mode. On the other hand, their harmonics will 
typically be dominated by Yq components that do not suf- 
fer from cancellation. Therefore, Rc will be higher for the 
harmonic of the (2,0) mode. This conclusion is independent 
of the the angle between the line-of-sight and the pulsation 
axis. 



1972; we 



In Table Al, the value of Rc for 2F4 is much larger 
than, e.g., that for 2F1 and 2F2. While we should caution 
that 2F4 could be blended with other combinations, this 
is exactly what one would expect from the example given 
above, given our identifications of Fl and F2 with I = 1 and 
F4 with 1 = 2 (Paper II). Note that one would not expect, 
e.g., F4+F1 to have a value of Rc larger than those for 
combinations of two 1 = 1 mode s, su ch as F2+F1, F3+F1 
and F5+F1 (as observed; Table Al). This is because the 
product Y§Y± has no Yq component and therefore suffers 
cancellation. 



Looking in detail at Rc for all combinations involving 
Fl and another 1 = 1 mode, e.g., 2F1, F2+F1, we find that 
the values are not consistent with being the same. There ap- 
pears to be a trend of increasing Rc with decreasing parent- 
mode amplitude. A similar trend appears to be present for 
the 3-mode combination frequencies. This is not expected, 
as generally the efficiency of mixing should be relatively in- 
dependent of the amplitudes of the input signals when these 
signals are weak. 

We should caution that above we have been neglecting 
the possible effect of different values of m. For instance, the 
harmonic of a f = 1, m / parent mode is described by 
a Y$ m distribution only; it would have a low value of /Jo- 
in our short observations, we cannot resolve the rotation- 
induced frequency splitting among different m components 
that has be en see n previously for some of the modes (Klein- 



et al. 



1998). Thus, the power at both the harmonic 
and its parent mode may result from a number of combi- 
nations, which makes our conclusions uncertain. They can 
likely be verified, however, using existing, much longer time 
series from WET. 

We now turn from these general considerations, valid for 
any local nonlinear mixing scheme, to a comparison of our 
data with numerical results and analytic calculations based 
on the convective-driving mechanism. As mentioned, in this 
picture flux perturbations appearing at the photosphere are 
diminished and delayed relative to those that enter the con- 
vection zone, by an amount that depends on the depth of 
the convection zone (and thus rc). The depth of the convec- 
tion zone, however, depends sensitively on the photospheric 
temperature. The latter varies due to the pulsations, and, 
therefore, the depth of the convection will vary as well. The 
associated variation in reduction and delay leads to a dis- 
tortion of the pulsation signal at the photosphere, which 
translates int o the p ower seen at the combination frequen- 
cies (Brickhill F992bl Wu [L997| ). 

Brickhill ( |l992b| ) found that R c will be ^6 for the first 
harmonic of an (l,m) = (1,0) mode and <; 16 for that of 
a (2,0) mode, the exact values depending on the frequen- 
cies, to and the angle between the line of sight and the 
pulsation axis. He found sum frequencies share similar val- 
ues of Rc with the harmonics^], while difference frequen- 
cies from the same pair of parent modes have lo wer v alues. 
These results were confirmed analytically by Wu ( Jl997j ), who 
found that for combinations with relatively high frequency 
((ui' ±3 tc 3> 1), -Rc should be approximately constant and 
A<E>c approximately zero, while towards lower frequencies 
{oj 3 tc — 1), Rc should be lower and A$c should shift 
towards —90° . 

The observed values of Rc and A"l>c compare, on av- 
erage, favourably to the expectations. The t ypica l values 
of Rc are reproduced with tc — 200 s (Wu 



1997), which 



is consistent with the values of tc inferred above from the 
longest-period o vers table mode and the velocity-light com- 
parison (Section 6.1). 



Little evidence is found for the expected trends of either 
Rc or A«l>c with frequency of the combinations. In part, this 
might be due to our inability to resolve components of dif- 



3 Brickhill (L992b) discusses but does not correct for the number 
of permutations n;j. 
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ferent m. Again, it would be good to verify these conclusions 
with longer time series, such as those already available from 
WET. 

As mentioned, a trend of Rc with parent-mode am- 
plitude is not expe cted in any simple non-linear mixing the- 
ory (e.g., Wu 1997). Nevertheless, the numerical calculations 



presented by Brickhill (1992b) show an amplitude depen- 



dence of Rc for 3-mode combinations that is similar to that 
observed (and with similar numerical values). At present, 
we do not understand why this dependence arises in the 
simulations; unfortunately, Brickhill does not expand on it. 



7 CONCLUSIONS 

We have applied the tool of time-resolved spectrophotome- 
try to the study of pulsating white dwarfs. In this paper, we 
showed that at our very high signal-to-noise ratios we can 
detect the small line-of-sight velocity variations associated 
with the g-mode pulsation s, at amplitudes similar to those 
expected (Robinson et al. 1982). Significant velocity varia- 
tion is seen for five out of the six modes we believed were 
real. For all five modes, the phase difference A3>v between 
velocity and light is between and 90°, i.e., velocity maxi- 
mum occurs later than flux maximum, but by less than the 
quarter cycle expected for adiabatic pulsations. The ampli- 
tude ratios Rv are similar for the five modes. The largest 
value occurs for mode F4 (at 776 s), which is consistent with 
mode F4 having 1 = 2 and the others £ = 1, as was inferred 
from the wavelength dependence of the pulsation amplitudes 
in Paper II. 

No significant line-of-sight velocity variation is seen for 
any of the five combination frequencies, even though two 
of these have light amplitudes comparable to that of the 
fourth strongest real mode. This provides independent con- 
firmation that the combination frequencies do not reflect 
physical pulsation. Rather, they likely result from mixing 
of mode power by a non-linear transformation in the outer 
layers of the star. The phase differences A$c between com- 
bination frequencies and parent modes are all close to zero, 
reflecting the sharp maxima and shallow minima in the light 
curve. The amplitude ratios Rc show significant differences, 
but are generally within a factor two of each other. The ex- 
ception is the four times larger value for the harmonic of 
mode F4. This difference is again consistent with the differ- 
ences in I between mode F4 and the others. 

Models for the interaction between pulsation and con- 
vection can broadly reproduce the observed values of Ry and 
A"I>v for the real modes, as well as tho se of R c and_A$c for 
the combination frequencies (Brickhill [l992b|; Wu Eg97]) . The 



models are internally consistent, in that the inferred ther- 
mal properties of the convection zone are such that modes 
at the observed periods are expected to be overstable. In 
detail, there are some problems, with A<E>v showing a trend 
with frequency opposite of that expected, and the mixing 
strength, reflected in Rc, not showing the expected corre- 
lation with frequency, but rather a possible dependency on 
parent-mode amplitude. Both these need confirmation; for 
the former, further time-resolved spectroscopy is needed, but 
for the latter, excellent data already exist. We also hope that 
our results will lead to further theoretical efforts. 

In summary, time-resolved spectroscopy offers insight in 



the white dwarf interiors at three levels. By measuring the 
effects of temperature variation at different wavelengths, the 
temperature structure of the atmosphere can be calibrated. 
Using the velocity changes and the combination frequen- 
cies, one can constrain the convection zone below the photo- 
sphere. And finally, with the determination of the spherical 
degree and hence the more secure identifications of period- 
icities with eigen-modes, one can use asteroseismology to 
probe the interior with greater confidence. 
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APPENDIX A: FURTHER PERIODICITIES 

The light curve we obtained is not described well by vari- 
ations at the eleven frequencies we identified in Section ^ 
In particular, there are a quite a few more peaks present in 
the Fourier transform at combination frequencies. We have 
continued the identification process until no peaks with am- 
plitudes larger than 0.3% remained. This gave us 29 frequen- 
cies^ At this point, the Fourier spectrum of the residuals 
became too noisy to identify with any confidence more peaks 
in the main region of power (1-4 mHz). In the cleaner region 
at higher frequencies, however, there were still four peaks co- 
incident with frequencies expected for 3-mode combinations. 
Since these might be of particular interest for modeling the 
pulsations, we included these in a final fit. The frequencies, 
ampli tud es and phases for all 33 periodicities are listed in 
Table |7~ 



Al 



Al i. We also fitted the veloc- 



, and indicated in Fig. 
ities to this set of 33 periodicities. The results are indicated 
calculated values of Rv 



Alb 



in Fig. 

modes are listed in Table |A1 



and A$v for the real 



We conclude with a few remarks. First, we note that 
possible significant periodicity is also present for 2F3 (Al = 
0.21 ± 0.05%, $l = -3 ± 14°, R c = 9.7 ± 2.4, A$ c = 
-83±14°) and F5+F2 (0.12±0.05%, -117±26°, 4.0±1.9, 
62 ± 26°), while no significant signal is seen at 2F5 (0.05 ± 
0.05%, Rc — 17 ± 17). One might have expected power at a 



number of combination frequencies not present in Table Al 
such as F3— F2 and F7±F1. These may indeed be present, 
but their frequencies are so close to other combinations that 
it is impossible to include them in the fit. Which of the 
possible combinations one uses is a matter of choice. We 
have been biased by theory in preferring those for which the 
parent modes have larger amplitudes, or, if similar, those for 
which the fitted value of A<E>c was closest to 0. 

Second, the periodicity F9, which is very close to Fl in 
period, is rather uncertain. The modulation could be fitted 
equally well with a variable amplitude for Fl. If F9 is real, it 
should have I ^ 1 given its proximity to Fl, which has I = 1. 
If it has 1 = 2 and a substantial m = component, one might 
expect a large harmonic, like for F4. This harmonic might 
in part account for the peak now described by F9+F1; if we 
fit 2F9 instead (we cannot fit both jointly), we find Al = 
0.26 ± 0.05%, $l = 109 ± 12°, which corresponds to R c = 



U. U.U. u. u. 

I I I I I ++ + + + + + ++ + + ++ + 



u 1 1 u/ w ii swin/ ii w/ 




2000 4000 
F (/zHz) 



6000 



Figure Al. Fourier transform of the light curve (a) and line-of- 
sight velocity curve (b), with 33 identified peaks indicated and 
labeled. The pluses indicate the amplitudes derived from fitting 
cosine waves (plus a constant for the light curve and a cubic 
polynomial for the line-of-sight velocity curve). The bottom curve 
in each panel shows the Fourier transform of the residuals, offset 
from zero for clarity. 



36±17, similar to that of 2F4, but A$ c = -160±20°, much 
further removed from zero than for any other combination 
frequency. 

Finally, while for most combination frequencies we do 
not find significant velocity variation, there are two excep- 
tions: F4-F1 at 2900 s and F6-F1 at 1447 s. These two 
periodicities appear harmonically related to each other, and 
we believe the associated velocity signals reflect either in- 
strumental or atmospheric effects rather than physical vari- 
ations. We briefly wondered about possible orbital modula- 
tion, but observations at higher resolution did not confirm 
these periodicities. 



^ The fitted amplitude is not always larger than 0.3% for com- 
bination frequencies, as for these the frequency is fixed at the 
sum or difference of the frequencies of the parent modes, which 
may not lead to a complete removal of the peak in the Fourier 
transform. 
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Table Al. Pulsation Frequencies in G 29-38 
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Notes. For the real modes, Ry = (Ay 


/2tt/)/Al in 


units of 


lOkmrad" 1 %" 


1 and A<J>v 


= *v - *L- 


For the 2- 



mode combination frequencies, / = f ± f 3 , with for difference frequencies i and j chosen such that / is positive, 



Rc 



In 



A i A 3 



with riij the number of possible permutations, and A<3?c = ^t, 3 ~ (^L-'-^l)' Similarly, for 3 



mode combination frequencies, / = /' + p + f k , R c = A L +j+fc ln^ h A\ j A\^A'l and A* c = ^ 3+k - (<I>^ + $J +$£). 
The uncertainties in the parameters for the light curve should be seen as indicating relative uncertainties only, as 
the light curve is not described well by the modes fitted and the deviations from the fit are definitely not distributed 
normally. The estimates given here were obtained by multiplying the observational errors by such a factor that 
-^red = 1' ^ or ^ e ve l° cl ty curve, the estimates of the uncertainties should be reliable. The unc ertainty distributions 
on amplitude and phase are not normal, but those on Acos<l> and Asin$ are (see Section 5.1). 
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